Non-Markovian Transmission through Two Quantum Dots Connected by a Continuum 
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We consider a transport setup containing a double-dot connected by a continuum. Via an exact solution of the 
time-dependent Schrodinger equation, we demonstrate a highly non-Markovian quantum-coherence-mediated 
transport through this dot-continuum-dot (DCD) system, which is in contrast with the common premise since 
in typical case a quantum particle does not reenter the system of interest once it irreversibly decayed into 
a continuum (such as the spontaneous emission of a photon). We also find that this DCD system supports 
an unusual steady state with unequal source and drain currents, owing to electrons irreversibly entering the 
continuum and floating there. 
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Coupling of a finite-size system to a reservoir with con- 
tinuum energy spectrum is relevant to many fundamental 
issues of physics, e.g., the emergence of classicality from 
a quantum world [[]]], and the origin of irreversibility in 
statistical physics. Simple example of such type of cou- 
pling can be a small system, with discrete states, coupled 
to a continuum in terms of quantum tunneling. 

In most cases, coupling with a continuum would result 
in an irreversible decay of the discrete state, with the par- 
ticle never coming back, such as the spontaneous emission 
of photon from an atom. However, in a recent study [||], 
where a double-dot connected by a continuum was con- 
sidered, it was found that the second dot can affect the de- 
cay from the first dot via quantum interference, and finally 
lead to a formation of stationary bound state embedded in 
the continuum. In subsequent related studies [|J Q], also 
concerning this dot-continuum-dot (DCD) system, it was 
demonstrated that the electron's motion through the con- 
tinuum is very unusual, say, revealing an undetectable fea- 
ture in the transfer process, and the information lost into 
the continuum can be retrieved without generating more 
disturbance. 

In the present work we consider inserting the DCD sys- 
tem into a transport configuration, as schematically shown 
by Fig. 1, which allows for an investigation for the behav- 
ior of continuous current through the DCD setup. More- 
over, beyond the above DCD studies [§-0], where a wide- 
band-limit (WBL) model was assumed for the central con- 
tinuum, we will generalize our study to a finite-bandwidth 
Lorentzian spectrum (FBLS) for both the leads and the 
continuum. As a consequence, the transient transport pro- 
cess is highly non-Markovian, much stronger than in the 
WBL [||]. Particularly, associated with the peculiar DCD 
setup, we will discuss and uncover other features in rela- 
tion with the coherence-mediated transfer and electron ac- 
cumulation (floating) in the continuum, basing both on an 
exact treatment for the non-Markovian transport. The is- 
sue of non-Markovian transport through quantum dots has 
received considerable attention in the past years [pl-O], 



where the non-Markovian transient dynamics are usually 
manifested in the current noise spectrum or in terms of 
full counting statistics. Moreover, exploiting the nature of 
non-Markovianity in broader contexts has been a subject 
of great interest in recent years, which stimulated novel 
means such as information-theoretic approach to recog- 
nize that non-Markovian dynamics is a reversed flow of 
information from the environment back to the open sys- 
tem [14 15]]. Also, experimental progress on revealing the 
non-Markovian effects was reported very recently | 
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FIG. 1: Schematic setup for the transport through double dots 
which are connected by a central continuum. 

Before solving the problem as shown schematically in 
Fig. 1, we first elaborate the method to be used in this 
work. As examined in Ref. [ [Tl| ] for the simplest problem 
of single-dot transport, either the usual Born-Markovian 
or non-Markovian master equation is not good in gen- 
eral for finite-band reservoir/lead, compared to the ex- 
act solution of this simplest model including also the re- 
cently derived exact master equation approach [pTft. In 
our present work, unlike using the Bloch-type rate equa- 
tion derived in Ref. Jl8| ] under the conditions of WBL and 
large bias voltage, we will apply an exact Schrodinger- 
equation-based single electron wavefunction (SEWF) ap- 
proach [|l9[ f^]. The SEWF approach is originated from 
an observation that, for noninteracting system (transport), 
the many-particle state of the whole system (i.e., the cen- 
tral system plus the two leads) can be constructed by 
the Slater determinant of all the single electron wave- 
functions. Note that, doing this, the Pauli principle is 
taken into account exactly. In practice, the single electron 



2 



wavefunction is solved directly from the time-dependent 
Schrodinger equation. Moreover, the initial single elec- 
tron states are summed over all the occupied states in the 
leads which are defined by the Fermi surfaces. In this 
way, the corresponding bias voltage is desirably intro- 
duced. Knowing the wavefunction of the entire system, 
of course, the current (or other quantities of interest) can 
be calculated by the matrix element of the corresponding 
operator and the total wave function. 

While in wide band limit, it is indeed possible to re- 
cover the master equation formalism from this SEWF ap- 
proach for large bias voltage Jl9[-pT|], but unfortunately so 
far similar conversion does not succeed for small volt- 
age despite the perfect numerical agreement with other 
exact approach (e.g., Ref. [|l7||). However, favorably, in 
practical manipulation the SEWF approach allows us to 
calculate the (transient) current directly from the single 
electron wavefunction, not needing to construct firstly a 
master equation. This is the most important advantage 
of the SEWF approach applying to noninteracting trans- 
ports. As we will see in the following, in practice the 
SEWF approach looks very like a generalization of the 
scattering-matrix approach of the Landauer-Buttiker for- 
mulation. That is, we first solve the single-electron state 
evolution and calculate its current. Then, we apply the 
well-known prescription to integrate the total current be- 
tween the two chemical potentials of the leads (the voltage 
window). So in certain sense the present approach can be 
regarded as a generalization of the Landauer-Buttiker for- 
malism, say, from the stationary to a transient version. 

Now we turn to solving the DCD system shown in Fig. 
1 by applying the SEWF approach. For simplicity, we 
assume that each dot has only one single spinless level. 
For noninteracting dot, where the spin-up and spin-down 
electrons transport independently, this model simply cor- 
responds to a single level between the bias window. In 
the presence of Coulomb interaction (inside the dot), this 
model is also reasonable for the case of a large Zeeman 
splitting so that within the (small) bias window there is 
only one energy level with definite spin (determined by 
the magnetic field orientation). In this case no Coulomb 
correlation inside the dot is relevant to the transport, since 
the dot level between the bias window can be occupied at 
most by only one electron. The whole setup is described 
by the following Hamiltonian 

H = y^ Eiajai + Eia\a\ + ^ E m a} m a m 

l m 
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Here, a\ {2) (a 1(2) ), a\ [r) (a l{r) ) and (a TO ) are, re- 
spectively, the creation (annihilation) operators of the left 
(right) dot, left (right) lead and the central continuum. The 



first four terms are the free Hamiltonians of the individual 
parts, while the last four terms describe their couplings, 
i.e., between the dots and the leads and the central con- 
tinuum, with the respective coupling amplitudes of fi; jT - jm 
and fL, ■ 



In this work, following Ref. (19[, we apply the SEWF 
approach to study the transient transport behavior. Let the 
initial state of the system correspond to filling the left and 
right reservoirs at zero temperature with electrons up to 
the Fermi energies /.il and p R , respectively. For nonin- 
teracting transport, the problem can be solved exactly for 
any values of the bias voltage, pi — fj, R . Indeed, the total 
wave function for the noninteracting electrons, ^(t)) = 
exp(— iHt)\^f(Q)), can be written at all times as a product 
of single-electron wave functions, \^{t)) = Yiil^ii^))- 
Note that, as well known, this product state in "second 
quantization" form is exactly equivalent to the construc- 
tion of a Slater determinant. Here \ipi(t)) describes a sin- 
gle electron initially occupying the state aj|0) in the left 
or right lead, with energy Ej = Ei n , We can express in 
general its state at time t as a superposition of all possible 
occupations in the whole system as follows 

\Mt)) = ^& i (i)a||0)+& 1 (i)a t 1 |0)+^6 m (0at n |0> 
+& 2 (t)4|0)+^& r (i)<4|0), (2) 



where |0) denotes the vacuum without any electron, and 
a ]|0) = 1) 2, Z,r,m) are the single-electron states with 
occupation amplitudes bj(t). Obviously, the assumed ini- 
tial condition is characterized by 6;(0) = %, while other 
amplitudes at the beginning are zero. Note that all the am- 
plitudes in Eq. (||) depend on I. However, for brevity, we 
did not express this dependence explicitly. 

For the single electron initially in state al\0), the 
probability of finding it in one of the leads or in 
the central continuum is therefore P^r m) fCO = 
Si( r ,m) \ b l(r,m)(t)\ 2 , and the average current reads 
I L{R,M),l{t) = cd t [P L(RiM)t i(t)}. Consider /ii > fj, R . 
At zero temperature, the total current in each lead or into 
the continuum can be expressed accordingly as 



Il(r,m) (*) 



(3) 



where is the density of states of the left lead. More 
explicit result for the single electron current Iur m) fCO 
will be presented in the following respective examples. 

Inserting the wavefunction, Eq. (^), into the time- 
dependent Schrodinger equation, we first obtain a set of 
coupled linear differential equations for the amplitudes. 
Then, applying the Laplace transformation, bj(uj) — 
Io° bj(t)e lujt dt, we convert the differential equations into 
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the following algebraic equations for bj (ui): 



at two equations for b\ (uS) and b 
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u~bi{io) - »6,(0) = EMui) + fi,&i(w)(4a) 

(w)(4b) 

i m 

w6 m (w) = E m b m (u>) + Cl m bi(Lj) + fj m 6 2 (w)(4c) 
wfo 2 (w) = & m b m {u) + E 2 b 2 {ui) + Y ^r6 r (w)(4d) 

m r 

u)b r (uj) = E r b r (uj) + f2 r &2(w).(4e) 

This set of equations can be solved as follows. First, 
find the solution of bi(uj), b r (uj) and b m (u>) in terms of 
&i(w) and b 2 (w), from the first, third and fifth equations. 
Then, substitute them into the second and forth equa- 
tions. In doing this, we will encounter the summation 
of Ei( r ,m) l^i(r,m) | 2 /[^ - E l(r . m) ], which can be further 
replaced by an integration by introducing the density of 
states, pu r ,m)' f° r me l ea ds and the central continuum. 
In the case of wide-band limit (WBL), we can regard 



Pl(r,m) 



and O 



l(r,m) 



as constants, accordingly denoted as 



Pl(r,m) and fl L ( RM ). We then have 



l^( r >™)l _ Xl(r,m) 



lj-Ei 

Z(r,m) 



(5) 
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where the constant coupling rates, ^l(r.m) = 
1-Kp L (R,M) \^l(r,m) I 2 , are introduced. 

In this work, of our greater interest is the more gen- 
eral case, say, finite-band Lorentzian spectrum (FBLS), 
for both the leads and the central continuum. This can be 
modeled as 
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/(r,m) 
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(0) 

L(R,M) 



where the constant £l£r R M \ and Kl(r.m) characterize, 
respectively, the spectral height and width. For the 
case ^l(r.m) much larger than the relevant energies, 
h-L(R,M) ^ w > the WBL is recovered. Otherwise, the 
coupling rates in Eq. (Q) are of energy dependence as fol- 
lows 
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L(R,M) 



Z(r,m) 



A 2 



(6) 
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ih-L(R,M) 
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where V 



(<>) 



io(0) 

L(R,M) ~ Z7r P L (K, M )\ il L(R,My 



Wide Band Limit. — In this case, from Eq. (Bh, we arrive 



ir L iT M 
u) - Ex + — - H — ) bi(u) 



= t— 62H (8a) 



- #2 + ~f- + "^J W 



(8b) 



Here we assumed Am = Qjvf and T^/ = f a/. While 
keeping i?i = —E 2 = e/2 7^ 0, we assume further 
= T R . = To and carry out an explicit analytic so- 
lution. From Eq. (^]), we solve for b~i(w) and b 2 {uS), and 
find that commonly each has three poles, Ei n and 0J\!2) = 
i(±A - 7) /2, where A = ^r 2 f - e 2 and 7 = T + r A/ . 
Then, straightforwardly, an inverse Laplace transforma- 
tion leads to 



b!(E in ,t) = 20i [Aie 

b2{Ei 7l ,t) = 2£Il£m 



-iE in t 



-«'/2ie 



/i2e 2t 



-iE in t 



/22e" 



(9a) 



(9b) 



In this solution, we introduced: /n = (2Ei n + e + 
%i)Z x Z 2 , /12 = [(A - ie)Z x + (A + ie)Z 2 ]/(2A), 
hi = ZxZ 2 , and f 22 = {Z 1 - Z 2 )/(2A), where Z 1(2) = 
e ±4t / a /[24 + i(7TA)]. 

With the knowledge of bi^(Ei n ,t), one is able to 
carry out the transient single-electron currents as fol- 
lows. As described previously above Eq. (||), from the 
time derivative of P L(RM) j = Ei(r,m) \ b l(r,m) | 2 , we can 
explicitly arrive to the following results: Ir,(-Ej„,t) = 
- eT L I &! (£ in , t ) 1 2 - 2ei7 L Im [e iE ^ b x (E m ,*)] for the left 
current; I R (E in ,t) = eT R \b 2 (E in , t)\ 2 for the right cur- 
rent; and 1 m (E in , t) = eT M \bi(E in , t) + b 2 (E in , t)\ 2 for 
the current flowing into the central continuum, from the 
both dots. Here we would like to stress that in steady 
state these currents satisfy the conservative relation 7j, = 
Ir + Iai, since the occupation probabilities in the two dots 
are stationary (unchanged) in the steady state. This con- 
servative relation can be checked numerically (in Figs. 2 
and 3) or proven analytically. The interesting feature is 
that, as we will see in the following, II 7^ I R even in 
the steady state! This is somehow a novel consequence 
of the continuum spectral property, but does not mean a 
breakdown of the electron number conservation. 

While some (complicated) analytic results are avail- 
able, in Fig. 2 we numerically display the transient single 
electron currents under the wide-band limit, representa- 
tively for an incident energy Ei n in alignment with the 
dot level. For simplicity, here and in the following numer- 
ical results, we consider two identical dots and symmet- 
ric coupling to the leads and the central continuum, with 
E x = E 2 and T L = T R = T M = f M = T a . 
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FIG. 2: Currents from an electron with incident energy in align- 
ment with the dot level (Ei n — Ei — E2 = 0), under wide- 
band limit for both the leads and the central continuum. 



Finite-Band Lorentzian Spectrum. — We consider ap- 
plying the finite-band Lorentzian spectrum (FBLS), 
Eqs. (g) and (Q), to both the leads and the central contin- 



uum. Again, for simplicity, we assume TV 

.(0) _ -p(O) 
M M 



(0) 



,(0) 



r£ } = f £? = T M , and A L = A R = A M = A. Straight- 



forwardly, solving Eq. for 61 (w) and b 2 (u>) gives 
2iQ in [e(u) + iA) + 2f3](cu + iA) 



b 2 (w) 



(w - Bj„)G(w) 
2iQ in AT M (u) + iA) 
(w - 



(10a) 
(10b) 



Here we introduced: 0,-„ = » / ^- s - 




A(r +r A/ )/2 + io;A,andG(w) = 4/3 2 - 
iA) 2 . For double dots with aligned levels, E\ — E 2 = 
e = 0, we find that 61 (u>) and 62 (w) have, commonly, 
five poles in the complex w-plane: Ei n , (±71 — zA)/2, 
and (±72 - iA) /2, where 71 = V2r A - A 2 and 72 = 
\J2TqA + 4r^/A — A 2 . For the purpose to express the 
transient current in a more compact form, we define in 
this context two types of inverse Laplace transformation: 



hit) = f&Me-**, and tyt) = 



with j = 1, 2. Then, after a short algebraic manipulations, 
the transient currents from a single electron with incident 
energy E in can be calculated by 

I R (E in ,t) = r AIm [b 2 (t)b* 2 (t)] (11a) 



-Tl (-E'in 1 1) 



-Im [r A6i(t)6I(t) + 2fi in &i(t)e ii3i » t l 



(lib) 



I M (E m ,t) = TMAlmlMt) + b 2 (t)][b* 1 (t) + &$(*)]} 

(11c) 

Here we omitted the energy in b 1 ^(Ei n ,t) and 
&*(2)(-^tru f° r tne sa k e of brevity. Integrating the re- 
sults of Eq. dTT| ) over Sj„ within a proper energy win- 
dow, transient currents for arbitrary bias voltage can be 
obtained. 

Before turning to the detailed behaviors of the currents, 
we first make a general remark on the transfer nature 
through the continuum. As we have mentioned in the in- 
troductory part, typically, a quantum particle cannot reen- 
ter the dots after it irreversibly decays into a continuum. 



To make this statement more specific, we may consider 
a single level (bound state) coupled to a continuum. In 
this case, the continuum can be a Fermi reservoir, but 
with the Fermi surface much lower than the bound-state 
level. Then, it is a common knowledge that an electron, 
initially occupying the bound state, will irreversibly decay 
into the reservoir along time, and will never come back in 
long-time limit. By contrast, for the case of DCD setup, 
the electron can transmit through the central continuum 
even in long time limit Below we explain that the ba- 
sic reason is owing to the coherent coupling of the two 
dots with the continuum. In this case, if we convert the 
above wavefunction approach/solution under into a rate 
equation, we will find that the (state) occupation transfer 
is mediated by a quantum coherence term, i.e., by the off- 
diagonal elements of the density matrix. This means that, 
if the "quantum coherence" is destroyed, the electron can- 
not propagate from one dot to another through the contin- 
uum. We demonstrate this feature, more transparently, as 
follows. In wide band limit, we can construct a Lindblad- 
type master equation for the coupling of the two dots to the 
continuum [||]: p = -i[H s , p]+T M V [a\+xa 2 ]p. Here p 
is the reduced density matrix of the double dots; the Lind- 
blad super-operator reads 2?[a]p = apa) — |{a^a, p}; and 
X = &m/£Im is introduced. We conclude that the co- 
herent superposition of a\ and a 2 in the Lindblad super- 
operator is of crucial significance to the charge trans- 
fer through the continuum, which actually provides such 
cross terms as a\pa\ to mediate the passage through the 
continuum. Otherwise, if we destroy the coherence by re- 
moving this type of cross terms and keep only terms as 
T>[ai]p and T>[a 2 ]p, we will find that the electron cannot 
transfer between the dots through the continuum. 
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FIG. 3: (a)-(c): Single electron currents for finite-band 
Lorentzian spectrum, with bandwidths of A = 0.2Fo, 0.5Fo and 
To, respectively. The electron's incident energy is aligned with 
the dot level which locates at the center of the Lorentzian. (d)- 
(f): Integrated currents over an energy window (— O.lFo, O.iro) 
around the dot level {E\ = E2 = 0). 

In Fig. 3 we show the result from the FBLS which, qual- 
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itatively, differs from the result under the WBL as shown 
in Fig. 2. That is, the current under WBL does not dis- 
play oscillations with time, while it oscillates in the case 
of FBLS. This oscillation behavior is indeed originated 
from an essential non-Markovian effect. The difference in 
Fig. 2 and Fig. 3(a)-(c) is out of a simple intuition: since 
the incident energy is in resonance with the dot level, we 
may expect the electron's motion to be affected mostly 
by the continuum states with similar energy of the inci- 
dent electron. However, by altering the bandwidth of the 
Lorentzian spectrum, the result in Fig. 3(a)-(c) indicates 
that the electron's motion is also influenced strongly by 
the "remote" continuum states with quite different ener- 
gies. 

The reason of having stronger non-Markovian behav- 
ior when coupled to a narrower band continuum may be 
understood in certain sense by imagining an extreme sit- 
uation: if the narrow band continuum is reduced to dis- 
crete levels in the same energy range, reversible quan- 
tum coherent oscillations will replace the oscillating de- 
cay behavior. There exists, however, an essential differ- 
ence between them. That is, no matter how narrow the 
continuum is, it will lead to an irreversible decay into the 
continuum. Of particular interest is that the electron is 
more easily coming back during the decay process, for a 
narrower bandwidth of the continuum. This understand- 
ing provides a dynamic insight for the nature of the non- 
Markovianity, which should be equivalent to the sophis- 
ticated information-theoretic notion by relating the non- 
Markovianity with a reversed flow of information from the 
environment back to the system of interest [|l4j |5j. 

We can further understand the stronger transient non- 
Markovian effect for a narrower band of continuum by 
means of the time-energy uncertainty relation. On rela- 
tively shorter timescale, the electron couples to a larger 
range of energies of the continuum. Then, for a wide 
band continuum, more "remote" energy states will involve 
into the (tunnel) coupling dynamics, resulting in a destruc- 
tive interference cancelation which destroys the "coming 
back" process (oscillating behavior). While on longer 
timescale (in steady state), the result is less sensitive to the 
bandwidth, as shown in Fig. 3(a)-(c). However, we find 
that the steady-state result under the WBL can be faster 
reached with a small increase of the bandwidth, while the 
short-time behaviors remain to be quite different. This 
feature demonstrates that the bandwidth effect hides in a 
more pronounced manner in the transient dynamics. In 
this context, we see that the transient single electron cur- 
rent captures more information than the conventional sta- 
tionary scattering approach. To make the result more rea- 
sonable, we further integrate the current over a small en- 
ergy window around the dot level and display the result 
in Fig. 3(d)-(f), where we find that the finite-bandwidth 
induced non-Markovian behavior is still evident. 

Finally we discuss the interesting steady-state behavior 
of the DCD system. The results in both Figs. 2 and 3 show 
that, even in steady state, the left and right currents are un- 
equal, i.e., Il > Ir- This feature distinguishes the central 
continuum, drastically, from any finite -size well, since in 



the latter case the steady-state left and right currents must 
equal to one another. The unequal left and right currents 
indicate a drastic transition when the finite-size central 
well approaches to a continuum limit. Actually, the es- 
sential physics is the development of irreversibility in this 
limit, which causes part of the transport electrons being ir- 
reversibly lost into the continuum and floating there, while 
other electrons supporting the transport current due to the 
coherent-coupling-mediated transfer between the dots. 

In this context one may expect a formation of chemi- 
cal potential in the continuum, owing to electron filling up 
there. Then, in steady state the net current into the contin- 
uum is zero. Of interest is that,in the central continuum, 
we did not introduce any inelastic scattering, which there- 
fore differs from what happens in the electrodes in trans- 
port, where the inelastic scattering processes (quickly) re- 
lax the electrode to a local equilibrium. For the central 
continuum in our DCD system, the absence of inelastic 
scattering implies that the electrons stayed in the contin- 
uum just float there, with no energy relaxation to form a 
local equilibrium with well-defined Fermi energy. We no- 
tice that, to account for phase breaking in quantum trans- 
port, Buttiker once proposed a type of dephasing model 
by introducing a virtual side reservoir [^2| j23|]. In that ap- 
proach, the transport electron first enters the side reservoir, 
experiences inelastic scattering in it, then returns back into 
the system. As a result, the quantum coherence of the 
electron is destroyed or partially destroyed. In Biittiker's 
dephasing model, the side reservoir is assigned an appro- 
priate chemical potential to guarantee a zero net current 
flowing into the side reservoir. On the contrary, in the 
DCD system, the nonzero net current flowing into the cen- 
tral continuum is a novel consequence of the continuous 
spectral property. 

Discussion and Summary. — Before summarizing the 
work, we would like to present a discussion on the non- 
Markovian nature of the SEWF approach. Indeed, con- 
ventionally, the non-Markovian feature is manifested in 
terms of a time-nonlocal master equation. The time- 
nonlocal feature, as is well known, reflects a memory 
effect. However, there is an alternative form of non- 
Markovian master equation, which is actually exact for 
some solvable models [O]. In that form, the non- 
Markovianity is manifested by a time-dependent decay 
rate, which depends on the entire past history thus keeps 
the full memory effect. Concerning the SEWF approach, 
it is built on a set of infinite coupled differential equations 
for the dot- and continuum-state amplitudes, see Eq. (Q). 
The basic solving procedure, in time domain, is integrat- 
ing first the continuum amplitudes, then substituting them 
into the equations for the dot-state amplitudes. In doing 
this, the dot-state evolution depends on the entire past his- 
tory of the continuum/reservoir, thus holds a memory ef- 
fect. This feature, in essence, corresponds to the time- 
dependent decay rate in the exact master equation formal- 
ism 0. 

While the non-Markovian behavior is rather obvious for 
the finite-bandwidth Lorentzian spectrum, it is of inter- 
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est to understand that, in the case of WBL, the transport 
under small bias voltage is in fact also non-Markovian, 
despite that the underlying non-Markovianity is usually 
quite weak, as indicated by It, in Fig. 2. Notice that, 
the conventional master equation is obtained by summing 
over all the single electron states initially occupied in the 
leads, which are actually defined by the Fermi surfaces. 
While for each individual electron, starting its evolution 
from a state in the lead, the transient dynamics is non- 
Markvian in essence as analyzed above, whether or not the 
master equation (obtained after the summation of the sin- 
gle electron states) is Markovian in WBL depends on the 
bias voltage. That is, under large bias voltage it is Marko- 
vian [18] while under small voltage, using the method in 
Ref. [17], it can be shown that the corresponding master 
equation is non-Markovian even in WBL. This latter re- 
sult is in agreement with the non-Markovian nature of the 
exact SEWF approach. 



To summarize, in this work we have presented an ex- 
act solution for the non-Markovian transport through a 
double-dot connected by a continuum. Unlike the intu- 
ition that a quantum particle does not reenter the dots after 
its irreversible decay into the continuum, we demonstrate 
a highly non-Markovian quantum-coherence-mediated 
transmission through the continuum, owing to coherent 
coupling of the two dots with the same continuum. The 
continuous spectral nature of the continuum also supports 
a steady state with unequal source and drain currents, as a 
result of electrons irreversibly entering the continuum and 
floating in it. 
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